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' , Abstract 

^. 
1) . 
,^ , We consider the production of intermediate-mass CP-even and CP-odd Higgs 

bosons in proton-proton and proton-anti-proton collisions. We extend the recently 

published results for the complete next-to-next-to-leading order calculation for a 

scalar Higgs boson to the pseudo-scalar case and present details of the calculation 

Q"^ \ that might be useful for similar future investigations. The result is based on an 

C^ ■ expansion in the limit of a heavy top quark mass and a subsequent matching to 

the expression obtained in the limit of infinite energy. For a Higgs boson mass 



^^ \ of 120 GeV the deviation from the infinite-top quark mass result is small. For 

300 GeV, however, the next-to-next-to-leading order corrections for a scalar Higgs 
boson exceed the effective-theory result by about 9% which increases to 22% in 
the pseudo-scalar case. Thus in this mass range the effect on the total cross section 
amounts to about 2% and 6%, respectively, which may be relevant in future precision 
studies. 



PACS numbers: 12.38.Bx 14.80.Bn 14.80.Cp 



Collider (LEP) 
the Tevatron [2]|1 



1 Introduction 

One of the most urgent problems in the modern particle physics is unveiling the origin 
of masses of elementary particles, which according to the Standard Model (SM) is closely 
related to the Higgs boson. A lower limit on the Higgs boson mass of about 114 GeV 
was set more than ten years ago by the experiments on the Large Electron-Positron 
and more recently mass values around 160 GeV were excluded by 
Due to rarity of the process, discovering the Higgs boson production 
requires subtle experimental methods and precise theoretical predictions. 

The Standard Model contains only one physical CP-even Higgs boson. However, many 
extensions of the SM, such as two-Higgs-doublet models or supersymmetric models predict 
also charged and CP-odd Higgs bosons. In this paper we consider the production of a 
pseudo-scalar Higgs boson in the form of an external current with a generic Yukawa 
coupling. We require that the coupling be proportional to the heavy quark mass, while 
the coefficient can be specified within any desired model. 

The dominance of the gluon- fusion process in the production of a scalar or a pseudo-scalar 
Higgs boson was established in the end of the 1970's [1H7]. Later, in the beginning of the 
1990's, several groups obtained the next-to-leading order (NLO) QCD corrections [814T0]. 
The latter appeared to be large, modifying the LO prediction by as much as 100%. Thus, 
the accurate prediction necessitated the next-to-next-to-leading order (NNLO) calcula- 
tion. 

The NLO results include the complete dependence on the partonic center-of-mass energy 
(s), the Higgs boson mass (M$) and the top quark mass (Mj) [inilll]. Assuming an 
infinitely heavy top quark, one considerably simplifies the calculation while introducing 
only a few per cent error. The NNLO corrections were first computed in this limit, for 
the scalar Higgs boson in Refs. [T2] - [T5] and for the pseudo-scalar in Refs. [T514T7] . More 
recently, the missing top mass-suppressed corrections to the scalar Higgs production were 
estimated and found small [T8H20] compared to other uncertainties. In these papers 
the mass dependence was recovered by interpolating between the expansion of the cross 
section near the heavy top limit and the leading asymptotics in the s — )■ oo limit, given 
in Ref. EOlEII. 



In this paper we apply the technique used in Ref. [19] and provide similar top mass- 
suppressed corrections to the pseudo-scalar Higgs boson production. We compute five 
expansion terms of the cross section in 1/M^, and match the partonic cross sections to 
the values of Ref. [22], derived in the s — )■ oo limit. 



Similarly to the conclusions of Refs. [TBH2U] . we find that the infinite-top quark mass 
approximation with factorized exact LO mass dependence receives relatively small cor- 
rections. We are not aware of any trivial explanations of this behaviour. Note that 
assuming no factorization of the exact LO quark mass dependence the 1/M^ result au- 



^We refer to [3] for critical comments on the Tevatron analysis. 



gumented with s — ?> oo behaviour deviates far from the infinite-top mass result, and only 
after including at least the l/M^ corrections the agreement of hadronic cross sections 
reaches the level of a few percent. 

For completeness let us also mention several results that improve upon the fixed-order 
QCD. Those include the soft-gluon resummation to next-to- next-to-leading [23] and next- 
to-next-to-next-to-leading |2lH26] logarithmic orders and the identification (and resum- 
mation) of certain vr^ terms [27] that significantly improves the convergence of the per- 
turbative series. Recent numerical predictions of Higgs boson production in gluon fusion 
both at the Tevatron and the LHC are summarized in Ref. [281430] . For reviews, see 
Refs. 



The remainder of the paper is organized as follows: in the next Section we introduce our 
notation and the basic formalism. After that, we describe the treatment of 75 for the case 
of the pseudo-scalar Higgs boson. In Section[3]we concentrate on the NLO prediction and 
compare our approximations to the exact result both at the partonic and the hadronic 
levels. Section |4] is devoted to the NNLO partonic corrections and Section |5] discusses 
the hadronic cross section. Results are presented for the LHC running at 14 TeV center- 
of-mass energies. Conclusions are presented in Section |6l In Appendix |A] we discuss the 
phase space master integrals expanded to the e order sufficient for a N^LO calculation. 
Technical details of the convolutions of various functions are given in Appendix [B] 



2 Preliminaries 

2.1 Notation and the LO result 

In the full theory, the scalar and the pseudo-scalar Higgs bosons couple to fermions via 
the following terms in the Yukawa Lagrange density: 

^y = -9V<^^<i'-9V<^^^i'<i'^ (1) 

where the dimensionless couplings g^'^ and g^'^ specify the coupling strength of the 
Higgs bosons to the heavy quark q. The cross section of the pseudo-scalar Higgs boson 
production is proportional to {gY'^Y'-, in the following discussion of this cross section we 
drop this constant for convenience. 

In the Standard Model g^'^ = 1, g^'"^ = 0, but in the MSSM, e.g., one has g^" ~ 1/ sin/3, 
g^j" ~ l/cos/9, g^" ~ l/tan/3, g^^' ~ tan/3 where tan/3 is the ratio of the Higgs 
field vacuum expectation values. Thus, the Higgs coupling to the top quark mass is 
suppressed for large tan /3 and our analysis is only valid for small values; for larger values 
also the contribution from bottom quarks has to be considered (see, e.g., the recent 
publications [331IM] ). 

The superscript "0" in Eq. ([T]) indicates bare quantities. Since we only consider QCD 



corrections, there are no counterterms for f°, if° and yl°, and the Higgs field vacuum 
expectation vahie is t" = 2~^/^G^ . We also introduce the variables 
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where s is the partonic center-of-naass energy, $ denotes either scalar {H) or pseudo-scalar 
[A) Higgs boson and Mt is the top quark pole mass. The divergences in the loop integrals 
are regularized using dimensional regularization with d = 4 — 2e space-time dimensions. 

The partonic cross section is commonly factorized as 
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where the factor Af^Q contains the constants and the complete non-trivial LO p- 
dependence. Such factorization provides an excellent agreement between the exact and 
the approximated results for the hadronic cross section at the NLO and the NNLO [T81420] . 
The well-known LO result is 
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In what follows, whenever we refer to the infinite-top quark mass result we assume the 
factorization of the exact LO mass dependence as given in Eq. ([3]). 

2.2 Optical theorem and asymptotic expansion 

Already at the NLO one has to consider real and virtual contributions which individually 
contain quadratic poles in e. In our approach we consider the forward scattering amplitude 
and use the optical theorem in order to derive the inclusive total cross section for the 
production of Higgs bosons (only the cuts that cross the Higgs boson line should be 
considered). At the NLO, the possible initial and final states are (7(7 — > $ + (0 or 1 gluon), 
qg ^ ^ + q, and gg — > $ + (7. At the NNLO, we in addition have reactions gg — )■ 
$ + {gg or qq), qg ^ ^ + qg, qq ^ ^ + {gg or qq),qq — )■ $ + gg, and gg' — )■ $ + gg'. Here g 
and g' stand for different massless quark flavourso Sample diagrams of the corresponding 



^It is understood that the ghosts are always considered together with gluons. 
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Figure 1: Sample forward scattering diagrams whose cuts correspond to the LO, NLO 
and NNLO corrections to gg^^ + gg,qg^^ + qg and gg -t- $ + qq. Dashed, curly and 
thick (thin) solid lines represent Higgs bosons, gluons and top (light) quarks, respectively. 

forward-scattering amplitudes are shown in Fig. [TJ 

The expressions are simplified by the forward-scattering kinematics implied by the optical 
theorem. The proper projectors applied to the external massless particles reduce each 
amplitude to a scalar expression depending only on p and x. Yet, the imaginary part of 
the double-scale four-loop integral with the exact dependence on both variables is still 
out of reach with the present methods, and we apply the asymptotic expansion [35] in 
p — > which corresponds to the limit M^ 2> M|, s. After the expansion, every four- loop 
integral factorizes into one-, two- or three-loop vacuum bubble with the single mass scale 
Mt, and the tree-level (for the virtual corrections), one- or two-loop box graph depending 
on X. 

The integrals with various denominator exponents that appear during the asymptotic 
expansion are reduced to a few master integrals using the integration-by-parts (IBP) 
relations, in which we treat cut lines as normal propagators. The two- and three- particle 
cuts are again re-introduced in the master integrals and evaluated separately. 

In the case of the virtual corrections the imaginary part of the Feynman diagrams arises 
solely from the factor 
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(a) (b) (c) 

Figure 2: Possible Feynman diagrams for the virtual corrections to the ^f^f-initiated par- 
tonic cross section. The shaded blobs represent massless one- or two-loop diagrams, the 
black dots mark the effective vertices that have no imaginary part. 



where in our case a = 1 or a = 2. Schematically, the occurring 
into three cases sketched for the (yf^f-initiated diagrams in Fig. 
(b) the calculation is straightforward, the loop integrals are 
of Eq. ([7]) expanded in e. If the grey blob develops an imag 
since it corresponds to a cut outside the Higgs line. The case 
integrals on the both sides of the Higgs boson propagator. Here 
corresponding product of a one-loop amplitude and a complex 
has to replace the factor (—1 -|- zO)^^ with (—1 -|- zO)^(— 1 — iO)" 



diagrams can be divided 
, 121 In the cases (a) and 
evaluated and the factor 
inary part, we discard it 
(c) has massless one-loop 

in order to reproduce the 
conjugate amplitude, one 
= 1. 



2.3 Treatment of 75 



The optical theorem simplifies the treatment of 7^ which appears in the coupling of the 
pseudo-scalar Higgs boson to quarks. We follow the prescription of Ref. [36] for the 
pseudo-scalar current renormalization, replacing 



75 -^ 



^e,,p.l^^YYY 



The square brackets denote total anti-symmetrization. The 24 terms on the right-hand 
side of Eq. (|8]) can be simplified via 



1 1 



24' ' " 4 ■ (^' 

Next, one should factor out the e tensor and compute the remaining amplitude in d 
dimensions. At the very end, when the expressions are finite in the limit e — )■ 0, one 
multiplies the result with the e tensor and applies the finite renormalization constant of 
Ref. [36]. 

Each forward-scattering amplitude contains two factors of 75. After the replacements of 
Eq. ([8]) we have the product of two e tensors that can be immediately re-written as the 
product of four metric tensors with antisymmetrized indices: 
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The right-hand side of this equation is defined in d dimensions and can be used during 
the calculation of the Feynman diagrams which eliminates explicit projectors. 



2.4 Alternative approach to virtual corrections 

As a cross-check, we evaluated the virtual corrections using two different approaches. 
First, we use the optical theorem as described above with a simple implementation of 75 
according to Ref. [36], treating the virtual and the real corrections on the same ground. 
The finite result is then obtained before applying the finite 75 renormalization constant. 

In the second method we consider the (pseudo-scalar) Higgs-gluon-gluon vertex diagrams 
and expand in the (formal) limit Mt ^ Mh- This is similar to the calculation of Ref. [37] 
for the scalar Higgs boson. The amplitudes are multiplied with projectors that couple 
to Lorentz indices of the gluons and the four additional indices that remain after the 
epsilon-tensor removal [38] . 

More explicitly, the pseudo-scalar Higgs-gluon-gluon amplitude has the form 

where a and a, /3 and h are the Lorentz and the colour indices of the incoming gluons. 
After the replacement Eq. ([8]) A'^g^A becomes 

^gg-^A,al3 ~ ^iiupcr^gg^A,al3 ■ \^^) 

Together with Eq. ( ITT1) it gives 

^gg^A,al3 = Wa^99-^A 9i I29a9p , (13) 

where gi and g2 are the incoming momenta of the gluons. Now the (Lorentz) scalar 
amplitude ^g^^A can be obtained from Ag^'^^^^ via the projector P^j^^. 



Aab paf3 Afj,vpcT,ab 

gg-^A — ^fiiypa^gg^A,al3 ' 

(d-2)id-3)iq,-q2y ^ ' 

In the actual calculations we use Eq. (|H]) in the initial diagram, then drop ef^,^pcr and 
obtain A'^g'^^^^^^. Using Eq. ( fT4l) and summing the diagrams, we then arrive at ^gg_s.A 
which has no open Lorentz indices. The virtual contribution to the total cross section is 
finally obtained by squaring the amplitude Eq. ( ITTl) and integrating over the phase space. 
Accounting for the averaging factors, we find: 

_ vr {d-3)\ab |2 ...X 



■^Note that the formulae presented in Ref. [38] apply to the coefficient function of the effective theory 
whereas here we invesigate the virtual corrections in the full theory. 



2.5 Software 

To evaluate Feynman diagrams we use the well-tested chain of computer algebra pro- 
grams developed for the scalar Higgs case [19]. The original diagrams, the corresponding 
asymptotic expansion prescriptions, and the master integrals are identical to the scalar 
Higgs case, with only trivial changes. Thus, here we only briefly outline the procedure. 

First, the diagrams are generated with QGRAF [39] supplemented by additional scripts that 
eliminate unnecessary graphs. Each diagram is then expanded in the limit M^ 3> s, M"^ 
using two independent programs: a combination of q2e and exp [IDISI], and a sepa- 
rate Perl program. This turns the original four-loop three-scale integrals into products 
of single-scale (one-, two-, and three-loop) vacuum bubbles and double-scale forward- 
scattering integrals (of one and two loops). The latter are reduced to master integrals 
with our own implementation of the Laporta algorithm |l2lll3]. Master integrals have 
been computed in Ref. [11]; in this calculation we extended them by one order in e which 
might be useful for the future calculations. 



2.6 Initial state singularities 

The renormalized sum of the real and the virtual diagrams is not finite in the limit e — > 0. 
The remaining poles originate from the coUinear divergences in the initial state: the initial 
gluon may split into a quark-antiquark pair, and the quark participate in the Higgs boson 
production. The corresponding contribution is determined by 

renormalized /^N fj r;i.barc <>, d i D /o, i-barel l^a\ 

-= K-[cTii, ®Pkj + Pik®cr^.^ J, (16) 

(5), bare 



^ij^'i>+X (^ 



TT 



where Pj • {x) is the splitting function that describes the probability of parton j to emit 
parton i with the fraction x of its initial energy. TZ is the renormalization operation, and 
ah ^^"^ is defined as prescribed by the definition of the splitting functions and consistently 
with the definition of the parton distribution functions (PDFs) which absorb the non- 
perturbative and non-singular features of the initial state. 

In our calculation, we use the common definition consistent with the MSTW08 PDF set, 
where ai ^^^ is the bare coupling in the effective theory with decoupled top quark and 
Uf = 5 massless quarks. The convolution of functions that enter into Pj • (x) and cr^j'^'^ is 
defined as 

[f®g]{x)= / dxidx2S{x -xiX2)f{xi)g{x2). (17) 

Appendix [B] contains the details of evaluation of these integrals. In this context, see also 
Refs. [T^[T5] and references therein. 



2.7 From partonic to hadronic cross sections 

In order to find the hadronic cross section one has to convolute the partonic cross section 
with the PDFs which can be written in the form 
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where the sum includes all distinct production channels, ij G {gg, qg, qq, qq, qq'}. The 
luminosity function dCij/dx contains the symmetry factors and the convolution of PDFs. 
For example, the quark-gluon luminosity is given bjQ 
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In Ref. [12] the luminosities dCij/dx have been discussed in some detail; for our purposes 
it is important to remember that dCgg/dx is practically zero for x < 10~^. 

We split the hadronic cross section according to 

(Jpp'_^$+xlsj = cr$ +c)cr^ +0a^ , (20) 

and add subscript "00" or "exact" to the quantities in Eq. ( 120|) when referring to the 
infinite-top quark mass result or the exact expression, respectively. A subscript "n" 
indicates that corrections through order 1/M" have been included. We leave out the 
subscript $ in case when the meaning is clear from the context. 

We use LO, NLO and NNLO PDFs by the MSTW2008 collaboration gl] in order to obtain 
the respective predictions for cTpp/^^^x in Eq- (I2D])- The choice of the PDFs determines 
the values of a,(Mz) = af\Mz): 

a^^{Mz) = 0.139384, ^^^^(M^) = 0.120176 , af^^O(M^) = 0.117068 . (21) 

The appropriate MS beta function then determines asifiR) that enters the formulae. 

3 NLO cross section 

3.1 Partonic cross section 

In this subsection we discuss the partonic NLO cross section as a function of x = M^/s. 
As mentioned in the previous section, we computed the partonic cross section in the limit 



"^This definition applies to pp collisions at the LHC; its modification for pp collisions at the Tcvatron is 
obvious. 



M^ S> s^M"^. By construction these approximations poorly converge for s > (2Mt)^, i.e. 
in the energy region where top quark pairs can be produced. In terms of the variable x, 
the threshold is given by Xthr = M|/(4M^^) which is a:thr = 0.12 for M$ = 120 GeV and 
a^thr = 0.75 for M$ = 300 GeV. For x < Xthr we do not expect our expansion in l/M^ to 
converge. 

To cure this problem we follow the suggestions discussed in the literature and match 
the heavy-top expansion to the leading term in the high-energy expansion obtained in 
Refs. [201 [21] and [22] for the scalar and pseudo-scalar Higgs boson, respectively. The 
matching procedure has been successfully applied to the scalar Higgs boson in Refs. [T8l - 
|20] . In the pseudo-scalar case the matching was done in Ref. [22] based on the leading 
order approximation in M^ ^ s, M\. In the following we discuss the effect of the power- 
suppressed terms and compare to the exact result as implemented in HIGLU 



Before we discuss our results in detail, let us describe our matching procedure. The 
expansion near the heavy top limit of the quantities Ajj (x) in Eq. (|3]) has the forrro 

A^f (x) = A,j- o(x) + p^^J,l{x) + p'A,,- 2(x) + . . . (22) 

and is expected to converge for x > Xthr- 

The a; — )■ limit of the NLO partonic cross section has the form 

A,,(x)"=°a, + 0(x), (23) 



with coefficients Cij computed in Refs. [20H22] . In order to combine these results we 
find some matching value Xm ~ a^thr and coefficient Dij, so that Aij{x) = Cij + DijX, 
X < Xm, and Ajj(x) = A°J''(x), x > Xm- As a matching condition, we require that 
^Ij^i^m) = Cij + DijXm and ^A^JP(xm) = Dij, i.e. we ensure that the transition between 
the a; — )■ and Mj ^ s, M^ approximations is smooth. The matching prescription for 
the channels with quarks in the initial state is less intelligent due to slower convergence 
below threshold and the oscillating functions. Here we set Xjn = 0.9a;thr and determine 
Dij from the requirement A'^^^{xm) = Cij + DijXm- 

In Fig. [3] we show the partonic cross sections for the pseudo-scalar Higgs boson production 
for Ma = 120 GeV (left column) and Ma = 300 GeV (right column). The approximations 
obtained in the limit of large top quark mass (dashed lines) show good convergence below 
Xthr and diverge, as expected, for small values of x. For the gg and qg channel the 
matching procedure leads to good approximations to the exact results for both values of 
Ma- The minor differences in the gg cannel only lead to small deviations at the hadronic 
level. On the other hand the deviation of the matched result to the exact one in the qg 
channel leads to a shift of about 15% in the hadronic contributions as discussed below. 
At the LHC with ^/s = 14 TeV the contribution of the qg channel to the NNLO part of 
the total cross section varies between about (—8)% and (—17)% for Higgs masses between 
110 GeV and 300 GeVJj Thus, an accuracy of 15% at NNLO induces an uncertainty of 



'^In what follows we omit the superscripts "$" and "(1)". 

^These numbers arc basically identical for a scalar and pseudo-scalar Higgs boson. 
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Figure 3: Partonic NLO cross sections for the gg (top), qg (middle) and qq (bottom) 
channel as functions of x for Ma = 120 GeV (left) and Ma = 300 GeV (right column). 
The expansion for p — )■ (dashed lines) is compared with the exact result (solid lines). 
Lines with longer dashes include higher order terms in p. The interpolation results (see 
text) are shown as dotted lines which (for gg and qg) demonstrate an irregular structure 
since several lines are plotted on top of each other. The threshold values Xthr are 0.12 and 
0.75, respectively. The subscript "qb" corresponds to the qq channel. 

at most 3% in the NNLO piece and is hence negligible in the sum of the LO, NLO and 
NNLO contribution. 

The situation is different for the qq channel. The partonic cross section has a charac- 
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+ 0.0002 
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exact 




15.4981 




16.6956 




16.1360 




19.1723 



Table 1: Contribution to the coefficient of 5(1 — x) in the normahzation of Eq. ([3]) from 
the various expansion terms p" = {M^/M^)"". The renormahzation scale has been set to 



teristic peak at x ~ a;thr which is reproduced neither by the x -^ nor by the p — t- 
approximation. Due to the small contribution of this production channel (below 1% at 
NLO) the accuracy of the infinite-top quark mass result is sufficient. 

It is worth mentioning that the quantities A)^ •• shown in Fig. |3] for the production of 
a pseudo-scalar Higgs boson are very similar to the corresponding scalar ones. This is 
both true for the shape and the numerical size of the corrections. For the gg channel it 



would not be possible to distinguish the curves for A 
a visible difference only occurs for x — )■ 0. 



A,(i), 
gg 



X 



and Af/'^ 



X, 



for X > 0.001; 



It is interesting to separately look at the ^-function contribution which is not shown in 
Fig. [3l In Tab. [T]we present the coefficient of 6{1 — x) in the normalization of Eq. ([3]) 
both for the scalar and pseudo-scalar case. One observes that the bulk of the contribution 
is given by the leading term whereas the power-suppressed terms add an additional part 
of only 8% and 17% even for M$ = 300 GeV in the scalar and the pseudo-scalar case, 
respectively. We anticipate that this behaviour is different at NNLO where the higher 
order terms provide more sizeable contributions as discussed in Section |H 

It is interesting to note that for M$ = 300 GeV the sum of the ffist five expansion terms 
in Tab. [1] reproduces the exact result to better than 1% in the scalar and 2% in the 
pseudo-scalar case. For the latter expansion, the error can be accounted for by doubling 
the contribution of n = 4 term. 



3.2 Hadronic cross section 

In this subsection we compare the hadronic cross section obtained from our approxi- 
mate expansions with the one obtained using the exact partonic results, implemented in 
HIGLU |l5] and integrated with the same PDFs. 



In Fig. m^a) we plot the exact LO cross section as a function of M$ for 200 GeV < M$ < 
400 GeV in order to demonstrate the threshold behaviour. Whereas there is a smooth 
transition in the scalar case, one observes a strong enhancement of the cross section for 
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Figure 4: (a) LO cross section for the production of a pseudo-scalar (top) and scalar 
(bottom) Higgs boson, (b) ratio of NLO and LO prediction for the total cross section. 
For these plots the exact results from [15] have been used. 

the pseudo-scalar Higgs boson. In Fig. ID^b) the purely NLO part of the cross section is 
shown where the LO contribution is divided out. This is actually the quantity which is 
approximated by our procedure. No strong enhancement is visible for Ma ^ 300 GeV 
which is the upper limit considered in this paper. Thus we expect that in the pseudo-scalar 
case the matching procedure works as well as for the scalar one. 

When discussing the quality of the matching procedure it is convenient to consider ratios 
of cross sections. We normalize either the exact or the matched expressions to the (un- 
matched) infinite-top quark mass result since these kind of ratios are also possible at 
NNLO. In the following plots the dashed lines represent the matched results including 
successively higher order in 1/Mt when going from short to long dahes. The solid line 
corresponds to the exact result. 

In Fig. [5](a) the NLO gg part is shown where the factorized expressions are used in the 
numerator. For Ma = 120 GeV the deviation between the various curves is below per 
cent level. For Ma = 300 GeV, however, one observes a deviation between the exact 
result and the infinite-top mass approximation of about 6%. The exact result and the 
matched result based on the leading order deviate by about 8% whereas there is perfect 
agreement with the matched result including 1/Mf terms. (The corresponding dashed 
curve is below the solid line.) 

An alternative comparison between the matched results on the one hand and the infinite- 
top mass and exact result on the other hand is shown in Fig. [5](b) for the gg channel where 
in the numerator of the considered ratio the fully expanded result is used and successively 
more terms in l/M^ are included. A large deviation of the matched result including only 
the leading Mt dependence (short dashes) from the exact result is observed. However, 
after including higher order mass corrections in the matching procedure good agreement 
up to Ma ~ 250 GeV is observed. For Higgs boson masses close to 300 GeV a significant 
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Figure 5: Ratio of the NLO parts of the hadronic cross sections for the gg (Fig. (a) 
and (b)), qg (Fig. (c)) and qq (Fig. (d)) channel. In (a) the exact LO result is factored 
out also in the matched result in the numerator whereas in (b) this option is abandoned. 

deviation of the matched and exact result is visible which can be traced back to the 
expansion of the LO result. Thus at NNLO we adopt the approximation of Fig. ^a) and 
factor the exact LO result out after including the l/Mf mass corrections. 

Figs, ^c) shows the result for the qg channel. The matched result agrees with the exact 
one with an accuracy of about 15% whereas the infinite-top quark mass result is about 
a factor two smaller than the prediction based on HIGLU. For ^/s = 14 TeV the qg part 
amounts between (-2)% (Ma = 110 GeV) and (-7)% (Ma = 300 GeV) of the total NLO 
part. Thus, if at NNLO a similar behaviour is observed it is important to incorporate the 
matched qg channels in precision predictions. 

As expected from the above discussion of the partonic cross section it is not possible to 
obtain a good approximation to the exact result for the qq channel which is shown in 
Fig. [5]^d). For lower pseudo-scalar Higgs boson masses there is also a significant deviation 
between the exact and the infinite-top mass result, for Ma = 300 GeV, however, they 
(accidentally, as can be seen in the lower right panel of Fig. [3]) agree perfectly well. 
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Figure 6: Partonic cross section for the gg channel for Ma = 120 GeV (top) and Ma = 
300 GeV (bottom). The plots on the left employ a logarithmic and on the right a linear 
x-scale. The notation is adopted from Fig. |3l 

4 Partonic NNLO cross section 

In this Section we discuss the partonic cross sections for the various channels at NNLO. 
We follow the matching procedure outlined in Section 13.11 for the NLO calculation, with 
minor modifications due to the different x — > limit. At the NNLOO 



Aij (x) ""= Eij In X 



0{1) 



(24) 



and Ref. [22] provides coefficients E^j. We again choose the matching value Xm ~ a^thr 
such that the transition between the x — )■ and M^ ^ s, M^ approximations is smooth: 
^ij''^(^rn) = Eij Inxm + Eij, ^A^JP(xm) = Eij/xm, and A'lJ^{xm) is given in Eq. ([22]). 

This prescription applies to the gg channel. For the other initial states we again use 
Xm = O.Qxthr and require that A'l^^^Xm) agree with the corresponding infinite-energy 
result. 

In Fig. [6]we show the partonic cross section for the numerically most important gg channel 
at Ma = 120 GeV and Ma = 300 GeV, with linear and logarithmic scales of the x axis. 



^We again omit the superscripts "4>" and "(2)'' 
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Figure 7: Partonic cross section for the qg channel for Ma = 120 GeV (left) and Ma = 
300 GeV (right panel). The notation is adopted from Fig. |3l 

The dashed lines represent the results for a'^J^ (cf. Eq. ( 12^ ): longer dashes indicate 
higher order p-expansion. As expected, these results diverge for x < Xthr- The results 
interpolated as described above are shown with dotted lines, shorter distances between 
the dots correspond to higher-order p-expansion in the matching procedure. 

For Ma = 120 GeV the partonic cross section has a maximum at x ~ 0.6 which is 
significantly higher than Xthr ~ 0.12 and thus it is nicely reproduced from a^J^^. The 
approximations start diverging at x ~ Xthr- As far as the matched results are concerned 
one observes stabilization starting from the one including the p^ terms. 

For Ma = 300 GeV we have Xthr ~ 0.75 which is close to the steep rise of the partonic cross 
section. The expansion results show good convergence properties down to x ~ 0.6 and 
start to diverge around x = 0.45, significantly below Xthr- As far as the matched results 
are concerned a similar behaviour as for Ma = 120 GeV is observed: The inclusion of 
higher order terms in p stabilizes the matching procedure, leading to firm NNLO results 
for the partonic cross section. 

Let us next discuss the g(7-initiated partonic cross section shown in Fig. [7] for Ma = 
120 GeV and Ma = 300 GeV. Again, very nice convergence happens at x > Xthr and the 
matching to s — )■ oo results is rather stable. The approximation is shown with dotted lines, 
the shorter distances between the dots denote higher-order terms in p. One observes that 
the latter are important for reliable results. Both for Ma = 120 GeV and Ma = 300 GeV 
there is a visible difference between the curves including only p'^ terms and the 0{p) 
result. Further corrections are relatively small. Judging from the behaviour at NLO it 
can be expected that the results for the partonic qg channel approximate the (unknown) 
exact result quite well with an uncertainty below 15%. 

The matching procedure was also applied to the quark-initiated processes qq, qq and qq'. 
We show the results in Fig. E] using the same notation as in Fig. El For the qq we expect 
a peculiar structure similar to the one at NLO (cf. Fig. E]) which is not reproduced by 
our procedure. However, since the overall contribution from this channel is small, we may 
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Figure 8: Partonic cross section for the qq ("qb"), qq and qq' ("qp") channels for Ma = 
120 GeV (left) and Ma = 300 GeV (right). The notation is adopted from Fig. El 

neglect the uncertainty of the matching procedure. 

The numerical contribution from the qq and qq' channels to the total hadronic cross section 
is similar to the qq one. For these initial states, however, we do not expect a complicated 
shape of the curves and thus the matched result through p^ terms can be trusted both 
for Ma = 120 GeV and Ma = 300 GeV. (The differences between the various matched 
results are again small at the hadronic level.) 

By analogy to Tab. (H in Tab. [2] we present the NNLO contribution to the 5- function part 
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pn 


Mh 


= 120 GeV 


Mh 


= 300 GeV 


Ma 


= 120 GeV 


Ma 


= 300 GeV 


n = 




83.5654 




94.0264 




95.9237 




105.0866 


n= 1 




+ 3.5586 




+ 24.0520 




+ 6.5340 




+ 44.7899 


n = 2 




+ 0.2190 




+ 9.1031 




+ 0.5330 




+ 22.5153 


n = 3 




+ 0.0164 




+ 4.2818 




+ 0.0495 




+ 13.0672 


n = 4 




— 




— 




+ 0.0050 




+ 8.1572 



Table 2: NNLO contribution to the coefficient of 5(1 — x) in the normahzation of Eq. ([3]) 
from the various expansion terms p" = (M^/M^)"'. The renormahzation scale has been 

set to M$. 



which is not contained in the plots discussed before. For Ma = 120 GeV the convergence 
is fast. The p terms still provide about 7%, however, already the p^ terms are below 1%. 
At Ma = 300 GeV a good approximation to the exact result requires terms through p^. 
Estimating the contribution of the p" terms for n > 5 as twice the p^ contribution, we 
expect the accuracy of about 8%. For Ma = 280 GeV this number goes down to 5%. 

Since the analytic results are quite lengthy we refrain from listing them here explicitly. 
They are available on request from the authors. 



5 Hadronic cross sections 

In this Section we present the hadronic cross section for the individual channels and 
compare in each case with the infinite-top quark mass approximation in order to quantify 
the accuracy of the latter. Although the results for the scalar Higgs boson have already 
been discussed in the hterature [T81 - I20] we present results both for the scalar and pseudo- 
scalar case for comparisons. For the numerical results we use the nominal LHC center-of- 
mass energy ^/s = 14 TeV; for ^/s = 7 TeV the results look qualitatively very similar. 

Let us again start with the (^^f-induced channel. In Fig. |9]the ratio 






(25) 



of the NNLO contribution to the total cross section is shown. In the numerator the 
matched partonic results are used and in the denominator is the bare infinite-top quark 
mass result. As before, the lines with longer dashes include higher expansion terms in p. 
Note that for the scalar Higgs boson terms through p^ are available whereas for the pseudo- 
scalar case we could compute even the p^ terms. Actually, in the practical calculation it 
turns out that the diagrams with pseudo-scalar couplings lead to significantly fewer terms 
at the intermediate steps. We believe that this is due to elimination of many terms with 
anti-symmetric properties of the e-tensors remaining after the 75-matrices. 
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In the left column the exact LO result is factored out both in the numerator and the 
denominatocl. One observes that the correction for the pseudo-scalar Higgs boson is bet- 
ween 2.5% for Ma = 120 GeV and about 20% for Ma = 300 GeV which is significantly 
larger than for the scalar Higgs boson, reaching at most 9% [T81 - I20] . Judging from the 
difference of the two consecutive curves, we observe in both cases good convergence of the 
p-expansion even for Ma = 300 GeV. This gives us some confidence that the corrections 
including the highest powers in p are approaching the unknown exact result. 

The right panels in Fig. |9] show the ratio (125|) with the numerator expanded in p without 
factoring out the LO mass dependence. The plot for the scalar Higgs boson reproduces 
the results in the literature [T8H20], demonstrating that after including higher orders in 
p the ratio becomes less dependent of M$. The curve including p^ terms deviates from 
unity at most by about 5%. The situation is different for the pseudo-scalar case. For low 
Higgs boson masses the corrections also converge against a few per cent, for higher mass 
values, however, they amount up to about 10% demonstrating that the infinite-top mass 
result accompanied by the matching procedure is not sufficient to approximate the exact 
result with 1% accuracy or better. 

Similarly to the NLO case shown in Fig. |5]we observe also in Fig. [9] that the convergence 
becomes poor for Higgs boson masses above 250 GeV if the exact LO mass dependence is 
not factored out. Thus, for our final prediction we use the approach of the upper left plot 
in Fig. [5] where the effect of higher order terms in p (beyond p^) can be safely neglected. 

The hadronic cross section for the quark-induced channels, qg, qq, qq and qq' are shown 
in Fig. [To] with the notation similar to Fig. [91 In all cases we observe decent convergence 
with higher p powers. From the NLO analyses of Section [3] we expect that the prediction 
for the qg channel agrees with the exact result to within about 15%. Note that there is 
a noteable correction from including the first p-suppressed term. This behaviour can be 
explained by the pattern in corresponding partonic cross section discussed in the previous 
Section (see Fig. [7]) . 

As far as the qq channel is concerned, the matched results coincide with the effective- 
theory prediction within roughly a factor of two to three for 120 GeV < Ma < 250 GeV 
which decreases to about 1.2 for Ma = 300 GeV. The results for the qq and qq' channels are 
shown for completeness in the bottom row of Fig. [TOl The deviation between the matched 
and the effective-theory result in the considered Higgs boson mass range is about 1.5 to 
2. For ^/s = 14 TeV the overall contribution from the qg, qq, qq and qq' channels to 
the NNLO corrections amounts to (-8)% to (-17)%, 0.1% to 0.2%, 0.08% and 0.3%, 
respectively. Thus, at the current level of accuracy it is certainly possible to use the 
infinite-top quark mass result for the predictions originating from the qq, qq and qq' 
channels. For Ma = 300 GeV the infinite-top quark mass approximation for the qg 
channel may lead to an uncertainty of about 2% in the NNLO prediction of the hadronic 
cross section. 



^The slight deviation from unity of the curve including only the p^ term is due to matching effects. 
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Figure 9: NNLO contribution to the hadronic cross section of the 5f(yf-channel normahzed 
to the infinite-top quark mass result. Top: scalar Higgs boson; bottom: pseudo-scalar 
Higgs boson. The left plots the exact LO result is factorized whereas in the right plots 
no factorization is used. 

It is instructive to look at the individual contributions of the NNLO pieces to the hadronic 
cross section. In Tab. |3]we show for M$ = 300 GeV the contributions from the 6 function, 
the plus distributions and the remainder. The infinite-top mass result is confronted with 
the approximations based on the matching that incorporates corrections of order p" {n = 
0, 1, . . .). One observes that both for the scalar and the pseudo-scalar Higgs boson the 
virtual corrections grow by almost a factor of two compared to the infinite-top mass 
values. In both cases the contribution from the plus distributions are quite small and 
the contribution from the 6 function amounts to about a quarter of the remainder. The 
latter dominate the corrections and amount to about 85% in the scalar and about 80% 
in the pseudo-scalar case. The difference between the best prediction and the infinite-top 
corrections amounts to 9% and 22%, respectively. 

For the M$ values below 300 GeV the convergence of all individual contributions signifi- 
cantly improves and our error estimate (given by the size of the last known term) is below 
2% and thus is negligible. 

Let us finally show results for the cross section of the numerically dominant (yffyf-induced 
contribution. In Fig. [TT] the LO, NLO and NNLO predictions for the total production 
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Figure 10: NNLO contribution to the hadronic cross section for the inclusive production 



of a pseudo-scalar Higgs boson induced by the ''qg'' 



'W, 



^qq" and "gg'" channels. 
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6{l-x) 


.• • •] + 


rest 


Mt —)■ oo 


0.363 


-0.066 


2.555 


"0" 


0.363 


-0.066 


2.544 


"1-0" 


0.093 


0.002 


0.054 


"2-1" 


0.035 


0.001 


0.034 


"3-2" 


0.017 


0.000 


0.034 


"3" 


0.508 


-0.063 


2.666 



A 


(5(1 -a;) 


• • •] + 


rest 


Mt —^ oo 


1.213 


-0.194 


7.728 


"0" 


1.213 


-0.194 


7.703 


"1-0" 


0.517 


0.013 


0.380 


"2-1" 


0.260 


0.006 


0.236 


"3-2" 


0.151 


0.003 


0.185 


"4-3" 


0.094 


0.002 


0.128 


"4" 


2.235 


-0.170 


8.633 



Table 3: Individual contributions to the quantity (5cr^^J;° (i.e. the gg induced part) in 



Eq. f l20|) for A/<j, = 300 GeV for a scalar (left) and pseudo-scalar Higgs boson (right). The 
first line corresponds to the infinite-top quark mass result and "z" represents the matched 
result including terms of order p\ The bottom line contains the best available prediction. 



cross section both for a scalar (left) and a pseudo-scalar Higgs boson (right) is plotted. 
The dotted and the dash-dotted lines belong to the exact LO and NLO prediction and the 
solid lines represent our best NNLO expression, i.e. the matched cross section containing 
the p^ and p^ corrections for the scalar and pseudo-scalar Higgs boson, respectively. For 
comparison we also show as the dashed line the result which is obtained if the infinite-top 
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Figure 11: Total cross section for the production of a scalar (left) and pseudo-scalar 
(right) Higgs boson to LO (dotted), NLO (dash-dotted) and NNLO (solid) precision. The 
dashed line corresponds to the prediction where for the NNLO corrections the heavy-top 
approximation has been used. 

quark mass approximation is used for the NNLO term. In the scalar case one barely 
sees a difference with the solid line whereas for the pseudo-scalar Higgs boson there is a 
deviation of about 6% for masses around 300 GeV. Thus it is important to replace the 
infinite-top mass results in this Higgs boson mass region by the matched results presented 
in this paper in future precision studies. 



6 Conclusions 

In this paper we consider the inclusive production of the scalar and the pseudo-scalar 
Higgs boson at hadron colliders. We compute the NNLO production cross section for 
Higgs bosons with masses below approximately 300 GeV, including top quark mass effects. 
The technique that we use extends the method that was recently applied to scalar Higgs 
boson production to the pseudo-scalar case. It relies on an asymptotic expansion in the 
limit of a large top quark mass and a subsequent matching to the zero- mass limit. This 
paper contains a detailed description of the method and several intermediate results which 
could be useful in the context of other calculations. In particular we provide all one- and 
two-loop four-point master integrals with an e expansion sufficient for a N^LO calculation. 

The main conclusions of our NNLO analysis are as follows. 

• At NLO our approach reproduces the exact result with high precision (see discussion 
in Section [3]). 

• Since at NLO the agreement between our approximation and the exact result is 
below 1% for the gg channel once expansion terms up to order 1/Mf are included 
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we are quite confident that the same is true at NNLO. Thus, from the practical 
point of view our calculation is equivalent to the exact one (at least for Higgs boson 
masses below threshold). 

• Our NNLO corrections deviate from the infinite-top quark mass approximation 
(with exact LO mass dependence factored out) by about 2% for low Higgs boson 
masses. For higher masses the deviation is larger. For M$ = 300 GeV it amounts 
to about 9% for the scalar Higgs boson and about 22% for the pseudo-scalar case. 

• This leads to the conclusion that up to an uncertainty of about 2% for the scalar 
case and about 6% for the pseudo-scalar case it is save to use the infinite-top quark 
mass approximation [T314T7J for the prediction of the total cross section. This is a 
non-trivial result which so far has no fundamental explanation. 

In case one aims for a better precision the infinite-top quark mass results have to 
be replaced by the results presented in this paper. 

• The accuracy of the NNLO part of the infinite-top quark mass result for M$ = 
300 GeV amounts to 6% (15%) for the scalar (pseudo-scalar) case at Tevatron. For 
the LHC with ^/s = 7 TeV the numbers are 8% and 20%. 

For Higgs boson masses above approximately two times the top quark mass the approx- 
imation used for the asymptotic expansion is not justified. In this limit one has to rely 
on the good agreement between the exact and the infinite-top mass result at NLO and 
assume that it extends to NNLO. 
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A Two-loop four-point master integrals 

The one- and two-loop four-point integrals have been studied for the first time in Ref. [1] 
(see Appendix B). Unfortunately, that reference contains a number of misprints. We in- 
dependently evaluated these integrals by a combination of soft expansion|j and differential 
equation methods. We furthermore extend the results of Ref. [H] by adding more terms 
in e which are required for a third-order calculation of the Higgs boson production cross 
section. The results can be downloaded in Mathematica format from the webpage ^6\ . 



''We acknowledge help with cross checks of the soft expansion by Robert Harlander and Kenial Ozeren. 
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The integrals published in Ref. [H] suffer from multiple typographical errors. We have 
also found a number of less trivial inaccuracies. In particular, the integral of Eq. (B.3) is 
off by a factor of 2, the right-hand side of Eq. (B.27) must read 



p2^-2. 



,3/ 



l°g(^) 4^3 _ 4C2log(^) + 4Li3(^) + i^^ -2\og\z) + 0{e) 



the right-hand side of Eq. (B.22) should begin with 



(26) 



V'z-\l - z)-^^ { Li2(l -z) + i^4^ + ... L (27) 



the two terms —4 log(l — z)^ — 8 of the 0{e^) contribution in Eq. (B.21) must be replaced 
with +41og(l - z)'^ + 8, etc. 



B Convolution of partonic cross section and splitting 
functions 

The singularities associated with the collinear radiation of quarks and gluons from the 
incoming partons are described by convolutions of the partonic cross sections (Jij{x) with 
the splitting functions Pjk{x) where a convolution of the two functions f\x) and g{x) is 
defined as 

[f ® g\{x) = / dxidx26{x-xiX2)f{xi)g{x2). (28) 

The functions that appear in the LO, NLO, and NNLO cross sections and LO and NLO 
splitting functions include combinations of HPLs of weight one, two, and three with the 



factors X, 1 — X, and 1 + x, and the generalized functions 5(1 — x) and 



ln''(l-x) 



It is 



l-x 

convenient to transform those functions to the Mellin space, where the convolutions turn 
into products of Mellin images: 

Mn [/(x)] = j x''-'fix)dx, (29) 

Mr,[[f^g]{x)] = M„[/(x)]M„[(7(x)]. (30) 

The Mellin transforms of various functions present in the NNLO Higgs boson production 
were studied in the literature |17]. We, however, decided to relate all required results to 
a limited set of Mellin transforms of HPLs with a certain maximum weight. 
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Let us consider the Mellin transforms of HPLs of weight zero and one: 

1 



Mn[l] = 


5 

n 


M„[Ho(a;)] = 


1 

9 ' 


M„[Hi(x)] = 


Si(n) 
n 


M„[H_i(x)] = 


i-ir 



(S_iH + ln2) + 



In 2 



(31) 



n n 

where Hi are the HPLs of weight 1 08]. The harmonic sums S...(n) are defined as follows: 



SW 



KAr^) 



5^/a(0Sg(z), Ui) 



4 = 1 



r", a>0 

(-1)^ z", a<0, 



(32) 



One may define the weight of those sums as the sum of the absolute values of their indices. 

Similarly to Eqs. f l3T]) . Mellin images of HPLs of higher weights contain harmonic sums 
of higher weights and various transcendental numbers originating from various HPLs 
evaluated at a; = 1. We choose the representations where each Mellin image contains only 
n as the argument of sums and integer powers of n in the denominators, and no products 
of sums (the corresponding algebra is discussed later). 

Some related Mellin transforms can be found by index shifting and integration by parts: 



X 



M 






dx" 



Mn+k [fix)] , 

x"/(a;)|J-(r2-l)M„„i[/(x)] 



(33) 
(34) 



In the latter relation, the boundary term has to be properly defined. In the limit x — )■ 
it vanishes since we always consider n higher than the order of the highest pole that 
f{x) may have at x = 0. If f{x) is some HPL, then its limit at x — )■ 1 may be either 
zero, a non-zero constant, or singular as In (1 — x). Such a singularity is present e.g. in 
Hi(x) = — ln(l — x), so that Mellin transform of ^Hi(x) = j^ does not exist. On the 
other hand, such functions can be regularized by turning them into a plus-distribution, 
and [jz^] has a well-defined Mellin transform, — Si(n — 1). We notice that this result 
can be obtained if we artificially drop the logarithmically divergent contributions from 
the boundary term in Eq. (134|) . In general, we may define the regularized derivative dx 
acting as follows: 



M„ 



dJix) =/?[/(x)]-(n-l)M„„i[/(x)] 



where 



R[gaix) ^'^(l -x)+ gb{x) ln''(l - x) + ... + goix)] = go{l) 
a,b,... > 0, gk{l) ^0 VA; > . 



(35) 



(36) 
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The extraction of divergent logarithms from HPLs is implemented e.g. in the function 
HPLLogExtract of the HPL.m package [I8] . 

Using this definition and the inverse Mellin transform, we may establish relations between 
the derivatives of HPLs and the common generalized functions: 



4Hi(x) 

4Hii(a;) 

4Hiii(x) 

4Hioi(x) 



5(1 -x), 





1 




1 


— X 

'ln(l 


+ 

-x)^ 




1 


- X 



1 

2 
6 



ln^(l-a;) 
1 — X 
1 



1 — X 



+ 5£iW^. etc. 

1 — X 



(37) 



Thus, it is not necessary to separately consider Mellin images of functions related to the 
derivatives of HPLs, such as -ji^, i_l ' ^^^■ 

At a given weight of sums, one may interpret the relations obtained by extending Eq. ( 13T|) 
as a system of equations, and complete it with the generalized derivatives of every line 
(excluding trivial relations). It is possible to solve this system for the monomials l/n*', 
S...(n)/n'^, and (— l)"S,,,(n)/?T,'^ in order to determine their inverse Mellin transforms. 
While writing this paper, we have learned about the work [19] that in particular mentions 
some very general algorithms to compute Mellin images of HPLs of arbitrary weights. 
Unfortunately, at this moment the program implementing those algorithms is not yet 
available to the public. 

If a Mellin image that we deal with is a linear combination of such known terms (with 
possibly shifted indices), its inverse transform is trivial to find with the help of the linearity 
of Mn[f{x)] and Eq. ( 133|) . However, to find the convolution of two functions we need to 
multiply the two Mellin images and the resulting terms may have a more complicated 
structure. Nevertheless, it is possible to transform the expressions to the required shape. 

First, the factorized denominators {n + a)~*(n + h)^K.. must be decomposed by partial 
fractioning. Second, the arguments of harmonic sums must be brought to agreement with 
the denominators using the definition Eq. (152]) : S^ g(n + 1) = S^ g(n) + fa{n + l)Sg(n + 1). 
Finally, the products of sums must be transformed into linear combinations of sums of 
higher orders. The corresponding algebra originates from the obvious identity 






A=l 



i=l 



(3^ 



which is applied recursively. For every product, it is then possible to arrive at the decom- 
position such as S_i,2('T-)S-i(r2) = 2S_i _i,2('^) + S_i_2,-i(?t.) — S_i _3(r7,) — '^2,2{n). 
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The above algorithm has been implemented as a Mathematica program. Since the HPL 
package [H] is not capable of computing integrals Eq. (!29|) with arbitrary n, the table 
of Mellin transforms (similar to Eq. ( I3T1) ) has to be pre-computed. The regular form of 
the relations Eq. ( I3T1) and the further results greatly simplifies the solution. We have 
been able to compile the table for the transforms of HPLs up to weight five, which should 
be sufficient to evaluate convolutions relevant to the NNNLO correction to the Higgs 
boson production. For that purpose one has to consider among other contributions the 
convolution of the NNLO partonic (^(^-induced cross section with the LO splitting function. 
The former contains contributions like Hqiq{x) / {1 — x) which can be written in terms of 
generalized derivatives of HPLs of weight 4. The LO splitting function has contributions 
from HPLs of weight 1 thus resulting in quantities of weight 5. Similary, the convolution of 
the NLO partonic cross section with the NLO splitting function involves HPLs of weight 2 
and 3, respectively, again leading to HPLs of weight 5. 
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